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O 

^vj , Abstract. This paper is concerned with a strongly degenerate convection- 

diffusion equation in one space dimension whose convective flux involves a 
r^ ' non-linear function of the total mass to one side of the given position. This 

equation can be understood as a model of aggregation of the individuals of a 
population with the solution representing their local density. The aggregation 
mechanism is balanced by a degenerate diffusion term accounting for dispersal. 
In the strongly degenerate case, solutions of the non-local problem are usually 
^^ ' discontinuous and need to be defined as weak solutions satisfying an entropy 

■ J* , condition. A finite difference scheme for the non-local problem is formulated 

^ , ■ and its convergence to the unique entropy solution is proved. The scheme 

emerges from taking divided differences of a monotone scheme for the local 
PDE for the primitive. Numerical examples illustrate the behaviour of entropy 
C^ , solutions of the non-local problem, in particular the aggregation phenomenon. 



oo 



^ 



1. Introduction 

> ; 

^D I 1.1. Scope. This paper is related to the initial value problem for a strongly degen- 

r^ ' erate convection-diffusion equation of the form 

Ut+U'U u{y,t)dyju{x,tyj = A{u),^, x e R, < f < T, (1.1) 

O; u{x,0)=uoix)>0, xeR, uo e (L^nL°°)(R) (1.2) 

for the density u = u{x, t) > 0, where A{u) is a diffusion function given by 

A{u):^l a(s)ds, where a(u) > for u e M. (1.3) 

^ ' The model (jl.ip , ()1.2p was studied as a model of aggregation by a series of authors 

including Alt JL^, Diaz, Nagai, and Shmarev [5], Nagai [H] and Nagai and Mimura 
[22l [23l [24], all of which assumed that a{u) =0 at most at isolated values of m. It 
is the purpose of this paper to study (11.11) , (jl.21) under the more general assump- 
tion that a{u) = on bounded u-intcrvals on which (jl.ip reduces to a first-order 
conservation law with non-local flux. We assume that 

A{u) — 7> oo as M — > oo. (1-4) 



C^ 
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2 BETANCOURT, BURGER, AND KARLSEN 

This implies that a(-) may vanish only on bounded subintervals of R. 

The key observation made in previous work [1] [2TJ [22l [231 El] is that if all 
coefficient functions are sufficiently smooth, and u(x, t) is an L^ solution of the 
problem (|l.ll) . (|1.2p . then the primitive (precisely, of u(-,i)) defined by 

v{x,t):= f ui^,t)d^, te(0,T], (1.5) 



is a solution of the local initial value problem 

v{x,0) = vo{x), a; G M, vo{x) 



G(0,T], 


(1.6) 


^ wo(OdC. 


(1.7) 



As a non-linear but local PDE, ()1.6|) is more amenable to well-posedness and nu- 
merical analysis. In this work we use that the transformation to the local equation 
(|1.6p is also possible in the strongly degenerate case, in which solutions of p.ip 
are usually discontinuous and need to be defined as weak solutions. To achieve 
uniqueness, an additional selection criterion is imposed, and the type of solutions 
sought are (Kruzkov-type) entropy solutions. The core, and essential novelty, of 
the paper is the formulation and convergence proof of a finite difference scheme 
for (|l.ip . (|1.2I) (in short, "w-scheme"). The scheme is based on a monotone differ- 
ence scheme for the initial value problem (11.61) . (jl.7p (in short, "w-scheme") in the 
strongly degenerate case, which in turn is a special case of the schemes formulated 
and analyzed by Evje and Karlsen [11] for the more general doubly degenerate equa- 
tion vt + ^{v)x = B{A{vx))x- The u-scheme is obtained by taking finite differences 
of the numerical solution values generated by the i)-scheme. 

The w-scheme is, in particular, monotonicity preserving, so the discrete approx- 
imations for V are always monotonically increasing when the initial datum vq is, 
and therefore the u-scheme produces non-negative solutions. Moreover, by modifi- 
cations of standard compactness and Lax- Wendroff- type arguments it is proved that 
the numerical approximations generated by the u-scheme converge to the unique 
entropy solution of (jl.ip . (|1.2I) . An appealing feature is that the primitive (II. 5p 
never needs to be calculated explicitly (except for the computation of v^). Nu- 
merical examples illustrate the behaviour of entropy solutions of (|1.1|) . (|1.2|) . and 
recorded error histories demonstrate the convergence of the v- and w-schemes. 

1.2. Assumptions. We assume that uq has compact support, and that there exists 
a constant A4 such that 

TY{uo)<M. (1.8) 

We also need that $ E C^(R), and that $ has exactly one maximum: 

3?;*eR: $'(«*)= 0, <^'{v) > ioi v < v* , <^'{v) < ior v > v* . (1.9) 

This assumption is introduced to facilitate some of the steps of our analysis; it is, 
however, not essential. In fact, in our convergence analysis of Section |4| we need to 
discuss the local behaviour of the numerical solution for v close to where it includes 
the value v* since that value is critical in the definition of the numerical flux. If 
we employ a function $ that has several separate extrema, then the locations of 
solution values including extrema are spatially well separated since the discrete 
analogue of Vx is bounded, and the techniques of Section |4| can be extended to that 
case in a straightforward manner. 
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1.3. Motivation and related work. Equation (|1.1|) . or some specific cases of it, 
were studied in a series of papers [U HI ETJ [22l [231 121] , in all of which it is assumed 
that a(-) vanishes at most at isolated values of its argument, so that it is always 
ensured that A'{u) > for m > 0. The interpretation of (jl.ip as a model of the 
aggregation of populations (e.g., of animals) advanced in those papers is also valid 
here and can be illustrated as follows. Assume that u(x,t) is the density of the 
population under study, and consider the equation 



ut+ { -k 



u{y,t)dy- u{y,t)dy u] == A(u)^^, fc > 0. (1.10) 

) Jx i / X 

Here, the convective term provides a mechanism that moves u{x, t) to the right 
(respectively, to the left) if 



u{y, t)dy < / u{y, t) dy (respectively, ...>...). 

) J X 

In other words, an animal will move to the right (respectively, left) if the total 
population to its right is larger (respectively, smaller) than to its left. Now assume 
that the initial population is finite and define 



Co := / UQ{x)dx. (1.11) 

It is then clear that (|1.10|) is an example of ()1.1|) if $'(u) = —k{2v — Cq), i.e., 

^(v) =-kv{v- Co) + const. (1.12) 

The aggregation mechanism is balanced by nonlinear diffusion described by the 
term A{u)xx, termed density-dependent dispersal in mathematical ecology. A typ- 
ical novel feature addressed by the present analysis is a "threshold effect", i.e. 
dispersal only sets on when the density u exceeds a critical value u^ > 0, i.e. 

I == ii u <Uc, Uc > 0, 
a[u) < 

I > ii u > Uc- 

More recently, spatially multi-dimensional aggregation equations of the form 

ut + "^ ■ {uW K * u) ^ AA{u) (1.13) 

have seen an enormous amount of interest, where the typical case treated in lit- 
erature is A = 0. Here, K denotes an interaction potential, and K * u denotes 
spatial convolution. For overviews we refer to [3l [HI [25] . The non-local and dif- 
fusive term account for long-range and short-range interactions, respectively, as is 
emphasized in [S]. The derivation of (I1.13P from microscopic interacting particle 
systems and related models, and for particular choices of K and A, is presented in 
[21 [4| [5l [HI [20] . Related models also include equations with fractional dissipation 
that cannot be cast in the form ()1.13|) . see e.g. [TFIFTT]. 

The essential research problem associated with (|1.13|) (or variants of this equa- 
tion) is the well-posedness of this equation together with bounded initial data 
u{x,0) = uo{x) for X € K'^, where d denotes the number of space dimensions. 
While the short-time existence of a unique smooth solution for smooth intial data 
is known in most situations, one wishes to determine criteria in terms of the func- 
tions K and A (or related diffusion terms), and possibly of uq, that either ensure 
that smooth solutions exist globally in time, or that compel that solutions of ()1.13p 



4 BETANCOURT, BURGER, AND KARLSEN 

will blow up in finite time. This problem is analyzed in ^2ll3ll4l[5l[6l[T6 l [T7 l [T8] 
(this list is far from being complete). 

The occurrence of blow-up was analyzed in terms of the properties of K for 
A = in [2j |3]; if iiT is a radial function, i.e., K = K{\x\)^ then blow-up occurs 
if the Osgood condition for the characteristic ODEs is violated, as occurs e.g. for 
K(x) = exp(— |x|), while for a C^ kernel this does not occur ^. Li and Rodrigo 
|16[ 117) consider this particular kernel and describe the circumstances under which 
blow-up occurs if the aggregation equation is equipped with fractional diffusion. 

The present equation (II. 1[) can be written as a one-dimensional version of (|1.13p 
only in special cases. However, and as was already pointed out in [22], (|1.10p can 
be written as 

ut + {uk*u)^^ A{u),j,^ (1.14) 

with the odd kernel K{x) — — fcsgn(a;). Equation ()1.14p becomes a one-dimensional 
example of ()1.13p if we observe that K * u = K' * u, where K' denotes the derivative 
oi K, if we choose the even kernel 

K{x)^-k\x\+C, (1.15) 

where C is a constant. We can write this as K{x) = —k{\x\) for K{r) = r — C. 
Suppose that one uses this kernel in the multi-dimensional equation (jl.131) . It is 
then straightforward to verify that in absence of dispersal {A = 0), the kernel (jl.l5|) 
satisfies the integral condition for blow-up in finite time, see [21. One result of our 
analysis is then that a strongly degenerate diffusion term A{u)xx, accounting for 
dispersal, is sufficient to prevent blow-up of solutions of (|1.13|) provided that the 
condition (|1.4[) is satisfied. In fact, in the context of aggregation models that are 
based either on (11.11) or on the more recently studied equation (|1.13p . the present 
work is the first that incorporates a strongly degenerate diffusion term, i.e. involves 
a function A{u) that is flat on a u-interval of positive length. So far, diffusion 
terms that have been considered in (jl.l[) degenerate at most at isolated u-values. 
Nagai and Mimura ^ studied the Cauchy problem for equation (|l.ip under the 
assumptions A{0) = 0, A'{u) > being an odd function. The initial function for 
the Cauchy problem in [22j is assumed to be bounded, non-negative and integrable. 
They prove existence and uniqueness of a bounded and continuous solution to the 
initial-value problem. In [23] the asymptotic behaviour of solutions to the same 
problem was studied for the specific choice 

A{u) = u", TO > 1. (1.16) 

It seems that the analysis of (I1.13P with degenerate diffusion has just started. Li 
and Zhang [18] study this equation in one space dimension for the diffusion function 
A{u) — u^/3, which degenerates at u = only. On the other hand, the numerical 
simulations presented herein show that under strongly degenerate diffusion, typical 
features of the aggregation phenomenon such as "clumped" solutions with very 
sharp edges [2S] appear. 

Finally, we comment that there is also considerable interest in the well-posedness 
and other properties of the local PDE ()1.6p (or variants of this equation) under the 
assumption that A{-) is an increasing but bounded function, an effect usually de- 
noted by saturating diffusion, cf. [7] and references cited in that paper. In this case. 
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which is exphcitly excluded by our assumption (|1.4|) , Vx (though not necessarily v it- 
self) becomes in general unbounded. It is at present unclear whether well-posedness 
of our nonlocal problem (|l.ip , (ll.2p can also be achieved under saturating diffusion. 

1.4. Outline of the paper. The remainder of this paper is organized as follows. 
In Section [2] we state the definition of an entropy solution of (|l.ip . (|1.2p . and point 
out that an entropy solution is also a weak solution. In Section [3. II we state jump 
conditions that can be derived from the definition of an entropy solution, and in 
Section 13.21 we prove the uniqueness of an entropy solution. 

Section |4] presents a convergence analysis for the u-scheme, which in part relies 
on standard compactness properties for the w-scheme. In Section HTTl the schemes 
are described. Section |4?2] contains a series of lemmas stating uniform estimates on 
the numerical approximations generated by the v- and the u-schemes, which allow 
to employ standard compactness arguments to deduce that both schemes converge. 

The final convergence result (Theorem 14. ip and its proof are presented in Sec- 
tion 231 This proof involves a discrete cell entropy satisfied by the w-scheme, which 
eventually permits to conclude that {m"_-^ /g} converges to an entropy solution as the 
discretization parameters tend to zero. This means, in particular, that an entropy 
solution exists. The mathematical model and the v- and u-scheme are illustrated 
by numerical examples presented in Section [5] 

2. Definition of an entropy solution 

Definition 2.1. A measurable, non-negative function u is an entropy solution of 
the initial value problem p.l|) . (II. 2|) if it satisfies the following conditions: 

(1) Wehaveue L'^inT)nL'^{nT)nL°°{0,T:BV{R))nC{0,T;L\R)), and 
A{u) €L'^iO,T;H\R)). 

(2) The initial condition (jl.2p is satisfied in the following sense: 

lim / \u{x, t) — uo{x)\ dx = 0. (2.1) 

(3) For all non-negative test functions Lp € C^{)1t), the entropy inequality 

JJut '^ (2.2) 

+ \A{u) - A{k)\tpxx\ dx dt > 

is satisfied, where v{x,t) is defined by (|1.5p and Ht := M x (0,r). 

Definition 2.2. A measurable function u is said to be a weak solution of the initial 
value problem (|l.ip . (|1.2p if it satisfies items (1) and (2) of Definition \2A\ and if 
the following equality is satisfied for all test functions (j) G C^{I1t)' 

|u(</.i + $'(w)0^) + A{u)(l)xx} dxdt = 0. (2.3) 

It is straightforward to check that an entropy solution of the initial value problem 
(|l.ip . (|1.2p is, in particular, a weak solution. 

Lemma 2.1. Assume that u is an entropy solution of the initial value problem 
(|l.ip - (|1.2l) (cf DeHnition \2. 1\) . Then u is a weak solution (cf. DeHnition \2.^] . 
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Proof. Choosing k > \\u\\ i^ao (^n^-^ in p.2p we obtain 

i-{u-k){(j)t + ^'iv)<f)^) - A{u)(j),:A dx dt > -k // u<^"{v)(j)dxdt 



'n 
or equivalently, 



|u(0f + <^'{v)(j)^) + A{u)(j)^^^ 



dxdt 



<fc// {(t>t + {<^'{v)4>)^\dxdt = 0. 

On the other hand, since we look for non-negative solutions, it suffices to set fc = 
in (|2.2p to deduce that we always have 

|u((/)t + $'(w)0x) + A{u)(l).^A dxdt > 0. (2.5) 

Combining (|2.4p and (|2.5p we see that u satisfies (|2.3p . D 

3. Jump conditions and uniqueness 

3.1. Rankine-Hugoniot condition and entropy jump condition. Assume 
that u is an entropy solution having a discontinuity at a point {xo,to) € Ht be- 
tween the approximate limits u~^ and u~ of u taken with respect to a; > xq and 
X < xq, respectively. Standard results from the theory of entropy solutions of 
strongly degenerate parabolic equations imply that such a discontinuity is possible 
only if yl(w) is flat for w e I(w"', u+) := [min{u~,u"'"},max{M~,u^}]. In that case, 
the propagation velocity of the jump is given by the Rankine-Hugoniot condition, 
which is derived by standard arguments from the weak formulation (j2.3p : 

s = ^^ ^ ^_ (^^'{v+)u+ - <^'{v-)u- - (A(w),)+ + {A{u),)-y (3.1) 

Here, (A(u)x)^ and {A{u)x)~ denote the approximate limits of A(u)x taken with 
respect to a; > a:;o and x < xo, respectively, and v^ and v~ denote the corresponding 
limits of v{x,t). However, since v{x,t) is continuous, we actually have v~^ = v~ , 
and the Rankine-Hugoniot condition p.ip reduces to 

..■.-(,fa.w)- '^'"''':-'^'-"- . (3.2) 

In addition, a discontinuity between two solution values needs to satisfy the jump 
entropy condition 

V/c G (minju", w^}, max{?i~, u"*"}) : 

^'{v+){u+-k)-{A{u)x)+ ^ ^ ^ <^'{v-){u--k)-{A{u)x)- (3-3) 

u — k ~ ~ u^ ~ k 

Taking into account v^ — v^ and (j3.2p . this reduces to 
Vfc e (min{M^, M^}, max{w^, w^}) : 

{A{u)x)+ ^ {A{u)x)+ - {A{u)x)- ^ {A{u)x)- (3-4) 

u^ — k ~ u+ — u^ ~ u^ — k 

In particular, if A{-) is flat on an open interval containing I(u",m+), then the 
double inequality p.4p is trivially satisfled. That $(u) is smooth, greatly simplifies 
the jump and entropy jump conditions. 
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3.2. Uniqueness of entropy solutions. The uniqueness of entropy solutions is an 
immediate consequence of a result proved in |14| (of. also |12) ) regarding continuous 
dependence of entropy solutions with respect to the flux function. More precisely, 
we have 

Theorem 3.1. There exists at most one entropy solution (according to Defini- 
tioi i2.1\] of the initial value problem (|1.1|) . (|1.2|) . Moreover, there exists a constant 
C ^ C (max |(f>'|) such that 

\\u{-,t) - u{-,t)\\^^^^^ < C\\uo - uoIIli(r), yt € (0,T], 

where u and u are entropy solutions with initial data Uq and Ug, respectively. 

Proof. Let u be an entropy solution of the problem 

Ut + {V{x,t)u)^^A{u),,, Vix,t):^'^' (J u{y,t)dyY 

with initial data u(0, x) = "^0(2;), and let u be an entropy solution of the problem 

Ut + {V{x,t)u),^^A{u)^^, V{x,t):^^'(j u{y,t)dy\. 

with initial data fl(0, x) — uq{x). According to |12|ll3j. keeping in mind that u and 
u are of bounded variation, i.e., u,u G L°°{0,T; BV{M.)), there exists a constant C 
such that 

II II f\ - I 

\\u{-,t) -u{-,t)\\^,^^^ < \\uo-uo\\l^r) + / \Vx{x,s) -Vx{x,s)\ds 

+ I \Vix,s)-V{x,s)\TV{u{-,s))ds 
Jo 

< ||wo-mo||li(R) +C' / \yx{x,s) -Vx{x,s)\ds. 
Jo 

Observe that 

/ |T4(a;, s) — t4(a;, s)| ds < max|<i>'| / \u{x,s) — u{x, s)\ds, 
Jo Jo 

so that by the Gronwall inequality we arrive at 

\\u{-,t) -u(-,i)|| 1 < exp(max|$'|i)||uo - wo||li(r)- 



D 



4. Convergence analysis of numerical schemes 



4.1. Preliminaries. We define the vectors U" := {u?j_i/2}jgz and V^ :— {w"}jgz, 
and discrctizc R by Xj :— jAx, j S Z, and the time interval [0,r] by i„ — 
nAt, n = 0,...,iV, At := T/N, N e N. We denote by u"+i/2 the ceU aver- 
age over Ij :— [xj,Xj^i] at time t„ and j g Z. We also define A :— At/ Ax and 
pL := At/ Ax"^ = A/ Ax and wherever convenient use the spatial difference operators 
A+0J := 0J+1 - 0j, A_(/)j := (j)j - (jjj-i, and 

A (/)j := A+A_(/)j = (j>j+i — 2(j)j + 4>j~i. 
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We assume that the initial datum uq is discretized via 



"°+i/2:=i / «o(OdC, jeZ. (4.1) 



Moreover, we define the operator S\x f^nd its inverse S^^ via 

5a.(C/";j) := Ax Y. "r+i/2, S-^liV'-^j) ■■- ^%7^- (4-2) 

/ — — oo 

Clearly, Sax and 5^;J; are the discrete analogues of the integral and differential 
operators that convert u(-,t„) into v{-,tn) and vice versa, respectively. Since we 
assume that ug is compactly supported, the sum in (14.21) is actually finite. 

The numerical scheme for the initial value problem (jl.ip , (11.21) can be compactly 
written as follows: 

C/"+i= [5X^oHo5ax]C/", n = 0,...,iV-l, (4.3) 

where the basic idea is to utilize a standard scheme of the form 

V''+^ ^n{V''), n^O,...,N-l (4.4) 

for approximate solutions of the local PDE (|1.6p , starting from the initial data 

^° := Ax Y, "?+i/2 = / ' "o(Ode, J e Z. 

/— — oo 

Clearly, if Cq is the total mass defined in (jl.lip . then we have that 

< w° < Co, w° < i;°+i for all J e Z. (4.5) 

Let us emphasize here that (j4.3p implies that 

C/" = [5^;^ o H o 5a J "f/° = [5a^ o H" o 5a J f/°. 

This means that for the actual computation of L/" from [/", the operators 5ax and 
'^Aa; need to be applied only once, and not for every time step. 

To derive properties of the scheme (|4.3p . we first analyze the scheme (|4.4I) . which 
is here given by the marching formula 

?;;+'=w;^-AA+[/i(w;_i,w;)-A(A_v;VAx)], j e Z, n = 0,l,2,..., (4.6) 

where A is subject to the CFL condition stated below, and 

/i(«;,z) :==$(0) + $+(w) + $_(z) (4.7) 

is the Engquist-Osher fiux [3], where we define the functions 



/ max{0,$'(s)}ds, $_(«):=/ min{0, $'(s)} ds 

JQ Jo 



We assume that At and Ax satisfy the CFL stability condition 

2Amax|$'(u)| +2umax|a(u)| < 1. (4.9) 

Note that the scheme for u can be written as 



u 



n+l 
J + 1/2 



u^^^,,^-\^+G^ + ^i^'A{u^^y.^), jeZ, n = 0,l,2,..., (4.10) 
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where we define 

^^ ■= Z^^+M-"-i.-") - ^ (^£ *VWd^ + £"' '^'_{s)ds^ . (4.11) 

For the ease of reference, we wih refer to (|4.6p - (l4.8p and (I4.10p . (|4.11l) as "f- 
scheme" and "u-scheme" , respectively. 

4.2. Uniform estimates on {«"} and {w"}. 

Lemma 4.1. Under the CFL condition (|4.9p . the v-scheme defined by (14.61) -i 
is monotone. 

Proof. We rewrite the scheme (|4.6|) as 

v;+'=n{v^_„v^,v^^+,)=:n], jeZ, n = 0,l,...,iV-l. 

Since a > 0, we then have 

(977" 

' -TAmin{0,$'(t;^±i)}+/^a(A±t-7/A2;) > 0, 



while the CFL condition (|4.9p implies that 

977" 

^ = 1 - A(max{0, $'(«;')} - min{0, $'(«;')}) - fiA+a{A^vpAa 



dv^ 



= I - \\^' {v'^)\ ~ ^lA+a{A_vy Ax) > 0. 



D 



As a monotone scheme, the scheme (|4.6p is total variation diminishing (TVD) 
and monotonicity preserving. Since (j4.6p represents an explicit three-point scheme, 
for a fixed discretization (Ax, At) we will always have 

Uj" = for j < -/C, w;' = Co for j > JC (4.12) 

for a sufficiently large constant /C > 0. Thus, we can state the following corollary. 

Corollary 4.1. // (|4.5P anrf t/ie CFL condition (|4.9p /lo/d, f/ien t/ie numerical 
solution {w"} produced by the v-scheme (|4.6p - (|4.8p satisfies 

0<f"<Co, f"<<+i for all j eZ, n = 1,...,A^. (4.13) 

yls a direct consequence, the numerical solution values V"' = {i'?}j"ez satisfy the 
(trivial) uniform total variation bound 

Tv(y") = ^|t;;-Vi-«"hCo. 

jSZ 

Lemma 4.2. T/ie numerical solution {w"} produced by the v-scheme (|4.6p - (|4.8p 
satisfies the L^ Lipschitz continuity in time property, i.e., there exists a constant Ci, 
which is independent of A := (Ax, Ai), such that 

m^r'-^j"! ^^1^- (4.14) 
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Proof. For j e Z, the quantity w" := v"^^ — v" satisfies 



n+3/2 n+1/2 , . r,/ „+l „+l 



-AA+[/.(.;.L+^^.r')-'^«-l'-")] 
AA+ [A{A^v'l+^ / Ax) - A{A^v'l/Ax)] . 



(4.15) 



We define 



I otherwise, 
and the quantities 

C;+'/':^[h{v^+\v]+l)-hiv^,v-+l)]e{v^+'~v^), (4.16) 

^n+i/2 _ [A{A+v';+^/Ax) - A(A+i;;/Ax)]0(A+z;;'+i - A+«;). 
Due to the monotonicity of {«"} (see (|4.13p ') we have 

c;+'/'>o, /?;+'/' >o, b;+'/'<o. (4.17) 

After some manipulations and using (|4.13p we obtain from (|4.15l) 

n+3/2 _ n+1/2 r-i ,^n+l/2 , nn+1/2 , / nn+1/2 p,n+l/2 



W, — W 



[1 - Ac;+'/' + ab;+^/' - A(i?;+i^/' + d;+'^')] 



Using the CFL condition we find 

l^n+3/2| < |^»+i/2| [1 _ A(c;+^/^ - b;+'^' + d;+1/' + d;+'^')] 

Summing this over j e Z, using (|4.17p and (|4.13p we obtain 

El n+3/2 1 ^ v^l n+1/2 1 

jGZ jGZ 

which imphes that 

El n+3/2 1 ^ Y^i 1/2 1 
iez jez 

From (|4.6p with n = we get 

jez jsz jez 

Using pr^ we arrive at (|i?Ti)) . D 



Lemma 4.3. T/ie numerical solution {«"} produced by the v-schemt 
satisfies the inequality \A^v^ / Ax\ < C3 with a constant C'3, which is independent 
of A. Equivalently, the solution {m"i_i/2} generated by the u-scheme (|4.10p . (|4.1ip 
satisfies the uniform L°° bound 

|w"+i/2| < Cs. for aU j G Z, n = 0, . . . , A^. (4.18) 
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Proof. It is sufficient to sliow tliat A{A+vJ / Ax) < C2 for a constant C2 tliat is 
independent of A. Taking into account ()4.12[) we get 

<\A{A+vpAx)~h{v^,vJ^,)\ 

= $(0)+ ^ A^{A{A+v'^/Ax)~h{v'^,v'^^,)) 



J 



k— — oo 
n+1 



k— — OG 



A 



$(0) 



<iEI<^' -^^1 + 1^(0)1 



Due to Lemma H:^ we see tliat \A{A+vpAx)\ < Ci if we choose C2 = Ci + |$(0)|. 
Taking into account (|1.4p concludes the proof. D 



Remark 4.1. Under all the previous assumptions we see that the v-scheme (14.61 
(|4.8p satisfies the following discrete cell entropy inequality 



,,n+l 



-k\- \v'; -k\+ AA+i/;_i/2 - Asgn(t;;'+i - k)A+A{A^v'; /Ax) < 0, 

where A; G K. anrf a numerical entropy flux is defined by 

H^+1/2 = ^;Vi/2 («", ^7+1 ■: k) := h{v^ V fc, «7+i V fc) - h{v^ A k, v^+, A k) . 

By a standard Lax- Wendroff-type argument we conclude that as A — >■ 0, the piece- 
wise constant functions assuming the value w" on Ij+1/2 x [imtn+i) converge to a 
limit function v that for all non-negative test functions (p £ C^{IIt) satisfies the 
entropy inequality 



V/ce 



||i; - k\ipt + sgn(t; - fc) [<P{v) - $(fc) - A{v^)]ip^^ 



dxdt > 0. 



However, since under the present assumptions v{-,t) is a monotone smooth func- 
tion, the entropy satisfaction property is not needed here. 

Lemma 4.4. The solution {w", j^ /j} generated by the u-scheme (j4.10p . (14. lip sat- 
isfies the following inequality, where the constant C4 is independent of A: 

TV{A{un) = Y.\A+A{u-_,/,)\ < C4. 



Proof. Using the marching formula (J4.6I) we can write 



\A+A{u-_,/,)\< 



1 



aI"^ 



n+l 



v';\ + \A+h{v^;_„v^)\ 



+ I {hiv-,v-) - h{v-_„v^)] e{v- - v-_,) \\A.v-\. 
Summing over j G Z yields 

jGZ jsz jez 

The right-hand side is uniformly bounded due to Lemma l42l and Corollary 14. II D 
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Lemma 14.41 does, in general, not permit to establish a uniform bound on the 
spatial total variation TV(f/") of the solution values {^"+1/2} generated by the 
w-scheme. This is possible only in the special case that A{-) is strictly increasing 
as a function of u, and a{u) vanishes at most at isolated values of u. 

We now prove that TV(C/") is nevertheless uniformly bounded, but by a bound 
that depends on the final time T. Our analysis will appeal to assumption (|1.9|) . 
From (|4.12p and (|4.13p we deduce that if {v^} is the numerical solution produced 
by the ti-scheme (|4.6p - (|4.8p . then at each time level there exists a unique index k 
such that w^ < w* < v^^^. The following lemma informs about the behavior of this 
index with each time iteration. 

Lemma 4.5. Assume that the data {v'^}j^z o-nd {«' }jez have been produced 
by the v-scheme ()4.6p - ()4.8|) starting from the m,onotone data {v'^Aj^j, under the 
CFL condition (|4.9p . Let k,k € Tj be the uniquely defined indices that satisfy 
Vk < V* < u^^.! and v^^^ < v* < w^iS respectively. Then k & {k - l,k,k + 1}. 

Proof. Since f^ < v* < v^^^ we analyze two cases: v^ < v* < v^^^ and v^^ < v* = 
v^^_^i. In the first, the monotonicity of the w-scheme and ()4.13p imply that 



Vk~^l S Wfc < W < Wfc+l S Wfc+2 ' 



such that either v'^+l < v* < v'^+\ or v'^+^ < v* < v'^+l, or v',^+1 < v* < v^+l, 
which means that fc = {fc — 1, fc, fc + 1}. In the second, we find that 

Wfc_l SVk<^ — "fe+1 ^ "fe+2 ' 

so either v^+l < v* ^ f^+\ or f^+i < v* < v^+l, or v'^+l < v* < v^+l We con- 
clude the proof by noting that w^^2 < ^* ^^ impossible due to the monotonicity of 
the w-scheme and (|4.13p . D 

The next lemma states the announced bound on TV(C/"). 



Lemma 4.6. Assume that the CFL condition (|4.9p is satisfied. Then there exist 
constants C5 and Cq, which are independent of A, such that the solution values 
[/" — {u^^,^,^}j^z satisfy the uniform total variation bound 

TV(C/") = ^|u;Vi/2-<-i/2| < (C5 + TV(t/°))cxp(C6r), n = l,...,N. 



Ji 



(4.19) 



Proof. From ()4.10p we obtain 



Let k be the index such that v^ < v* < v^^^^, and let us split Z into the subsets 



A 
B 
C 



:=^":={jeZ|j<fc-2}, 

:=6" :={jeZ|fc-2<j <A: + 2}, (4.20) 

:=C" :={.7 eZ|A: + 2<j}. 

Let w^ := A+u^_^ ,2 ^^^ a^' := A+yl(-u^'_^ ,2)6'(A+u^'_^ ,3). For j G yl, we obtain 
w';+^ = w}' - /iA_A2$(i;p + ^lA^[a'^w';). (4.21) 
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Using a Taylor expansion about w" we find that there exist numbers a" £ [w", wj+i] 
and /3" £ [w"_i, w"] such that 

Substituting this into (|4.2ip we obtain 

;A_(a>"(/3;')(A_.;f) 



+ A_$"(/?;) {A_v!;y + <i>"(/3;_i) (i^;-' - vj_^w:;_,Ax) 

+ 0(At) (u;^_i + w;.' + A+v'; + A^v]) . 
In an analogous way, we find for j £ C 

+ 0(At) (w'l + w;-Vi + A+w^' + A_ «;) . 

Now we deal with j £ B. For j = k — 1, using that v* is a maximum of $ and 
following analogous steps as before, we get 

<-i = <-i - m($KVi) - Hv*) + A_A2cf>«Li)) + MA'(aLi<_i) 
= <„i -/i($'(OKVi -i'*) + A_A2$(z;^0) +mA'K--i<-i) 
= <„i - n{{<P'{0 $'(^'*)) KVi - V*) + A.AH{v-^,)) 

+ /iA2(a;^i<„i) 
= u,I?_i [1 - A$'(«LO - 2^,al_,] + <_2 [Ma^2 + A$'(z;^2)] + ^^^a^ 

+ 0(Ai)(l + <_2 + <_i + A+<Li + A_^;^l)• 

For j = /c, using that $'(«*) = we compute 



+ /iA2(a>^) 
= < - ^[A+A^c&K) + A^A^c&K?)] + MA2(a>,") 

+ 0{At){l + <_i + < + <+i + A+v^^ + A^v-,). 
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For j = k + 1 and j = k + 2, the following steps are analogous to the prevfous cases. 
Using that $'(w*) = we obtain 



ti-\- 



I = <+i - A^[A+A2ci>(t,;Vi) + 3($K") - *(«*)) + *«) - $K"_i)] 



+ 0(At)(l + <+i + w^V2 + A+i^IVi + A_t;IVi)> 



w2t^ = w 



C'(At)(l 



"fe+2 + W^+3 + ^+Vk+2 



A^v 



k+2) 



Finally, summing over j we find that there exist constants Cq and C'l such that 



which implies that 



CV 



E K^' I ^ E l^°l exp(C6T) + ^ exp(ar), 
which proves (|4.19p . 



D 



The next lemma states L^ Holder continuity with respect to the variable t of the 
solution generated by (|4.10p . 

Lemma 4.7. The solution {w"i x/2} generated by the u-scheme (J4.10I) . (J4.11I) sat- 
isfies the following inequality, where the constant Cs is independent of A: 

y^hT+1/2 ~ "?+i/2|Aa^ <Csy/At{m - n) ior m > n, m,n e Nq- (4.22) 

Proof. We first establish weak Lipschitz continuity in the time variable. To this 
end, let (j){x) be a test function and 0^ := 4){jAx). Multiplying equation (|4.10l) by 
4>jAx, summing over n and j and applying a summation by parts, we get 



Ax^<^,(^;+i%-«j+i 



/2^ 



jez 



At^G7(</),-0,_i) 



+ ^Y.^<t>, - '/',-i) (^(";'+i/2) - ^(""-1/2)) 

Using Lemma 14.41 and the fact that is smooth we obtain 



AxY,h{^';tl/2~ ^1+1/2) <c\w\\At, 

iez 

where C is independent of A and 0. Consequently, for m > n the following weak 
continuity result holds: 



^^Y.^A<+i 



/2 "j + 1/2; 



j£2 



<C\\(p'\\At{m-n). 
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Since Ej :— u"^j^ /j — u'^^j^ has bounded variation on M, we arrive at the inequahty 
(|4.22p by proceeding as in [TOj Lemma 3.6]. D 

Now, following ideas of [13] we prove an L^ estimate for A'{-)x- 

Lemma 4.8. The solution {w", j^ /j} generated by the u-scheme (j4.10p . (14. lip sat- 
isfies the following inequality, where the constant Cg is independent of A: 

n=iiez \ / 

Proof. Multiplying (|4.10l) . by u"-,^,^Ax, summing the result over n = 0, . . . , A^ — 1 
and j G Z, and using summations by parts we get 

JV-l 

^EE(^-^K+i/2))(^--"+i/2) 

Af-1 . N-1 

n=0 jeZ n=0 j£Z 

A ^-1 

+T^EE(A-;v"i/2)^ 

n=0 jez 



where we used that 

("j + l/2 ^ "i+l/2)"j + l/2 = 2 V^J + 1/2) ~ ("j + 1/2) ~ (■"j+1/2 ~ "j + 1/2) 

In light of Lemma 14. 3[ we can also write 



since a(u) > 0. Using this observation, we find that 

N-l N-1 . 

i, E E(A-^KVi/2))' < At E Eg"(a-;Vi/2) + ^EKW)' 

n=0 jGZ n=0 jGZ jSZ 

A ^~1 

+f^EEK:i/2--"+i/2)'- 

(4.24) 

On the other hand, from (|4.10p and the inequality (a + 6)^ < 2a^ + 2b^ we obtain 

li^ltl/2 - ""+1/2)' < A2(A+G,")' + 2f,'[{A+A{u]^,/,))' + (A_A(.^+,/,))^). 
Multiplying the last inequality by Ax and summing the result over n and j yields 

^ E EKi% - "".1/2)' < ^ E E(A+Gj)^ 

n=0 jGZ n=0 jGZ 



JV-1 
n=0 jez 
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In what follows, we assume that the following strengthened CFL condition is satis- 
fied for a constant e > 0: 

CFL. := 2Aniax|$'(M)| +4uniaxa(w) < 1 - e. (4.25) 

The new CFL condition implies in particular that 

.2._,.At / At(l-e) 



ifi^Ax = 4^-— < 



Ax Ax a* 

and therefore 

N-l 



Ax 

n=0 jeZ 



EEi-ltlf 



(4.26) 



Summing (|T^ and (|i?^ yields 

< At E Eg"(a-;Vv2) + f^EK.1/2)' + ^ E E(a-g;vo' < ^, 

n=0 jeZ jeZ r!=0 JGZ 

where we used Lemma l476l the bound over G" and the fact that At = 0{Ax'^). D 

With the help of Lemma 14.81 we can prove 

Lemma 4.9. Under the assumptions of Lemma \4-8\ there exists a constant Cio 
which is independent of A such that 

^\A{u]\^/^) - A{u]^^/^)\^Ax < Cio(m - n)At for m > n. (4.27) 

jez 

Proof Using Lemma 1431 the fact that A'{u) > and (|4.10p we get 

E(^«+i/2)-^K+i/2))'Aa: 



< a* ^(A(u™ 1/2) - A(uJ+,/,)) («;•;,/, - z.;Vi/2) Ax =: ^ + S, 
where we define 

m—l 

A := -Ata*J2{A{uf_,,^,) ~ A{u-^,^,)) ^ A+G^, 

?n — 1 

B := Aa* 5](A(.™ ,/,) - A(.;Vi/2)) E ^^^K+i/^)- 
jez i=n 



(4.28) 
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Summing by parts, using the Cauchy-Schwarz inequality and Lemma l478l we obtain 



A = Ata*Y.J2 G]A^iM^^i/2) - ^K+1/2)) 

= At«* E E G\i^-M^T+i/2) - A-^K+1/2)) 

jez l=n 



m— 1 



<T«*EElG^I 



Ax 



Ax 



Ax 



m — 1 



+ ^a*J2J2 \G\\Ax^O{{m~n)At) 

j£Z l=n 

Proceeding in the same way for B yields 



13 = -Aa* J2\ [M^T+1/2) H^l+1/2) {A{<-i/2) - ^(-"-1/2))] 

m — 1 

/— n 

= -Aa* 5]|(A_A(..™ ,/,) - A-A(.7^,/,)) ^ A_^(^^,/, 

m— 1 

= -^«* E E (^-^«+i/2) • ^-M^Ui/2) - A_A(.^^,/,) . A_^(4^, 

CA-Afu" )\^ 
^^'+^^ Aa; = C'((m-ji)Ai). 
a; y 

Inserting into ()4.28p that yl, S = 0{{m — n)At) concludes the proof. 

Let us now denote by u the piecewise constant function 



/2 



D 



W-l 



'{x,t) ■- Y^ E^J"(^'*)"i+i/2' 
r!=0 jez 



(4.29) 



where Xjn denotes the characteristic function of Ij x [t„,i„+i), and let us denote 
by V its primitive. From the L°° bound (Lemma 14. 3p . the uniform bound on the 
total variation in space (Lemma 14.61) and the L^ Holder continuity in time result 
(Lemma 14. 7p we infer that there is a constant C such that 

\\u^\\L^(Ur) + h^\\L^in^)<C; h^(-,i)|B^(K)<C foralHe(0,T] (4.30) 

uniformly as Ax, At | 0, while Lemmas 14.81 and 14.91 implv that there are constants 
Cii and Ci2 independent of A such that 

(4.31) 



\A{u^{- + y, •)) - Aiu'^i; •))|L.(n,) < Cii(|y| + Ax), 
\Aiu^i;-+r))-A{u^i;-)) 



\L^(nT 



< Ci2\/'r + Ai. 
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4.3. Convergence to the entropy solution. 

Theorem 4.1. Assume that Ax and At satisfy the CFL condition (|4.25|) . and that 
uq is compactly supported and satisfies (II. 8|) . Then the piecewise constant solutions 
u generated by the u-scheme (J4.10I) . (|4.11l) converge in the strong topology of 
L^{I\-t) to an entropy solution of (|l.ip . (II. 2p (in the sense of Definition ] 2. 1]) . 



Proof Since w'^ G L°°{nT)r\L°°{0,T; BV{R))r\C^/'^{0,T;L\R)), we deduce from 
(|4.30p that there exists a sequence {AijigN with A, J, for i ^- oo and a function 
u€ L°^{nT)nL^{nT)nL°°{0,T;BV{R)) such that u^ ^ u a.e. ouHt. Moreover, 
in Hght of (|4.31l) we have A{u^) — > A{u) strongly on Li^^{IIt), and we have that 
A{u) S L^(0, T; i?^(K)). Lemma [4.71 ensures that u satisfies the initial condition 
(|2.H) . It remains to prove that u satisfies the entropy inequality ()2.2p . To this end, 
we show that the w-scheme satisfies a discrete entropy inequality, and then apply a 
standard Lax- Wendroff- type argument. From (|4.1ip we infer that 

where a^.^/j' /^i-i/2 *^ K-i'^^Fl satisfy 

<i>'^{a]^,/,) ^ 0{A+v-_,) r $;(s)ds, (4.32) 



^'_{l3-_,/,)^e{A+v-_,) $'_(s)d5. (4.33) 

Consequently, defining tlie function 

g-^,/,{u,v,w) ■.=uX<^'+{a]_,/,)+v[l - A(<i>V(a7+i/2) - ^'_{p-^,/,))] 
+ ^(-A<i>'_(/?;V3/2)) + K^i^) - ^Mv) + Aiw)), 
we can rewrite the scheme (|4.10p as 

"j + l/2 =^ »j + l/2(,"j-l/2''"j + l/2''"j+3/2J- 

Note that under the CFL condition, S". 1/2 is a monotone function of each of its 
arguments for all j G Z and n = 0, . . . , TV — 1, and that 






The quantity u'^lll/2 '■= '"?+i/2 ^ AA+G" + fiA'^A{u'^ -^^,2) satisfies for ah k G 



'^j+1/2 ^' ^ ^j + l/2("j-l/2'"j + l/2!'"j+3/2) ^j + l/2(^' ^'^) 

-Afc[A+$;(„;'_^/^) + A+ci>'_(/3;.Vi/2)], 

and since ^"1 j^/2 is a monotone function of each of its arguments, we get 
|-;;i% - k + A4A+$V(a;L,/2) + A+$'_(/3;-Vi/2)] | 

^ I^J + l/2("i-l/2''"i+l/2''"i+3/2) ^ Gj'+l/2ik,k,k)\ 
= ^7+1/2 ("i-l/2''"i+l/2''"j+3/2) '^ S^+l/2{k,k,k) 
- Sj+1/2 ("i-1/2) "^"+1/2. ""+3/2) A 5"+1/2(^, fc, k) 

< 5;+i/2 (""-1/2 V fc, w;Vi/2 V /C, l.;V3/2 V k) 

^ ^j"+l/2 ("i-1/2 ^ ^' "i+1/2 ^ ^! "j+3/2 '^ ^) • 
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Thus, defining 

G7(r,s,fc) := \r - m'+ia^^,/.^) + l-^ - fc|'i>'-(/3;Vi/2) 
-±-[\A{s)-A{k)\-\A{r)-A{k)\), 



we can write 



-n+l 



v:-,, - k + Afc[A+$V(a;L,/,) + A+$'_(/3;^^/,)] | 

On the other hand, 

\u]%, _ k + Afc[A+4>V(a^_,/,) + A+$'_(/3;^,/,)] I 

(4.35) 



Combining (|4.34l) and (I4.35p . we arrive at the "cell entropy inequality" 

l"j + l/2 '^l rj + 1/2 '^l +^^+*-^J l"j-l/2'^j + l/2''«j 

+ sgn(«;'^i/, - fc)Afc[A+<i>V(a^_i/,) + A+$'_(/3;.Y,/2)] < 0. 

We now basically establish convergence to a solution that satisfies (|2.2p by a 
Lax-WendrofF-type argument. Now, multplying the j-th inequality in (14. 36^ by 
Jj (p{x,tn) dx, where Ij := [xj,Xj+i] and (p is a suitable smooth, non-negative test 
function, and summing the results over j G Z, we obtain an inequality of the type 

El + E2 + E3 < 0, where we define 

w-i „ 

El := Y^ Xl(l"i+V2 ^ ^1 ^ Wj+1/2 - ^1) / ^{x,tn)dx, 

N-1 



E, := Afc Yl E^gi^K^i/2 - fc)A+(a>V(a^_i/2) + $'-(/3;Vi/2)) / Vix,tn) dx, 

n=0 jeZ -fj 

JV-1 
^3^=^EE^+^i("^l/2'"i+l/2'fc) / ^i^,tn)dx. 

„=o jez •'^3 

By a standard summation by parts and using that (p has compact support, we get 

Af-l 






<^(a;,i„+i) - ipj x.tn) 



n=0 jeZ J 

and £'3 = El + E"^, where 

N-l 

El := -M J2 E(l"^i/2 - ^l'J'V«i/2) + h;Vi/2 - fc|<5'-K+i/2 

n=0 jGZ 

/■ </j(a; + Ax, tn) - ip{x, t„) 

X / dx, 

Ji, Ax 

JV-l 

El:=-AtYY.\M^l+i/2)-^('')\^ 

n=0 j£Z 

f ip{x + Ax,tn) -2ip{x,tn) +'p{x - Ax,tn,x) 

X / — ^ dx. 

7/, Ax2 
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Clearly, we have that 

N-l 
n=0 

where 
1 k 

E2,n ■- ^ 






X 



N-l 

E" ^- ^ E E bsn{u]+l/, -k)~ sgn(^;Vi/2 - fc)] X 
n=0 jez 

J ij 

In the remainder of the proof, we will appeal to (|4.32p and (14.33^ , and assume 
that both $_ and $+ are smooth away from v* . Moreover, we know that for 
each n, the data {w"}jgz are monotone. Therefore we will utilize once again that 
for each fixed nG {0,...,7V— 1}, there exists k such that v"^ < v* < v^^-^ and using 
Lemma|13]we know that z;^+^ < v* < w^+j^ with ke{k-l,k,k + 1}. Thus, HA", 
B" and C" are the sets defined in (|4.20p . we may rewrite E2 „ as 

^2,n = ^2,yl" + ^2,i3" + ^2^0"- ' 

where the subindex denotes the summation over j from the sets yl", S" and C", 
respectively. For j £ A"" we note that 

A+<i>VK-i/2) + A+<i>'-(/3;+i/2) - A+ci>'(a^„,/2); 

using a Taylor expansion about v^ we can write 

A+^'4a]_,/2) = ^KVi/2 + u]_,/2)^"{v-) + 0{Ax'). 
Thus, we obtain 

EIa^^ E |sgn(u;Vi/2-^)(«"+i/2 + ""-i/2)'i>"(t'") / ^(x,i„)da; + 0(Ax). 
Since TV(f/") is uniformly bounded, this implies that 

EIa- ^kY, sgn(u;Vi/2 - ^)<+i/2'i'"(^") / -^(a^-irO dx + 0{Ax). (4.37) 
Analogously, we obtain 

4c" = k J2 sg^KVi/2 - ^)«;'+i/2*"(«,") / fi^^tn) dx + OiAx). (4.38) 

Moreover, since $ is smooth and we are inspecting the situation near the extremum 
(i.e. ^'{v*) = 0), we can conclude using Taylor expansions that 

A+$V(aJ_,/2) + A+$'_(/3;Vi/2) = OiAx) for j e S". 
Since ^B" is finite, S^^" = C'(Ax). Combining this with (|i371) and (|i351) we obtain 

4n-^E«g^("i+l/2-^)";Vl/2*"(«;) / V'(2:,in)dx + 0(Ax). 
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On the other hand, a summation by parts yields that 

Af-l 

1/2 



n=0 jGZ '- 

- A+$;(««_^/^) _A_<i>V(/3;Vi/2)} / ¥'(2^, in) da: 

+ [A+$V(a;'_i/2)+A+$;(/3;Vi/2)] / (<p(a;,t„+i)-(p(x,t„))da; 



Using arguments similar to those of the discussion of E^ and Lemma 14. 5[ we see 
that the expression in curled brackets is 0{Ax'^), and finally _E| = 0{Ax), so that 

N-l 

E, = kAt E E^gnK+i/2 - fc)«;Vi/2l'"(«") / Vix,t^) dx + OiAx). 

n=0 jeZ •'^3 

A treatment similar to that of E\ yields 



N-l 



^3^— AtEEI"".i/2-H*'(-") / ^^"^^"'y'^^"'^"^ dx+0(A.). 

Since (/? is smooth, we may state the inequality — Ei — E2 ~ El^ — E^ > Q as 

Uu^ - k\{ipt + '^'{v^)ip.^) - sgn(u'^ - k)u^k'^"{v^) 
+ \A{u^) ~ A{k)\ip^Adxdd > -CnAx 

with a constant Cn that is independent of A. Taking A — )■ we obtain that the 
limit function u satisfies the entropy inequality (j2.2p for almost all fc € M. To prove 
that (|2.2p is valid for all fc e M we may proceed according to Lemmas 4.3 and 4.4 
oi\m. D 



Remark 4.2. Theorem \4.1\ implies, of course, that an entropy solution exists. An 
inspection of the proofs of Lemmas \4-^ and \4-3\ reveals that the L°° bound for u is 
actually independent ofT. Thus, even though our analysis is limited to domains Wt 
with a finite final time T , we can say that entropy solutions of (jl.ll) . (II. 2p do not 
blow up in any finite time. 



5. Numerical examples 

The examples presented here illustrate the qualitative behavior of entropy solu- 
tions of the initial value problem (jl.ip , (jl.2p and the convergence properties of the 
numerical scheme. For the first purpose, we select a relatively fine discretization 
and present the corresponding numerical solutions as three-dimensional succesions 
of profiles at selected times or contour plots that should almost be free of numeri- 
cal artefacts, while the convergence properties of the scheme are demonstrated by 
error histories in some examples. For all numerical examples we specify Aa: and 
use Ai = At/Ax^ = 0.1, i.e.. At = O.lAa;^. 
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Figure 1. Example 1: Numerical approxim.ation of v (top) and 
corresponding approximation of u (bottom), obtained via (|4.3p 
with Aa; = 0.001. 



5.1. Example 1. In Example 1 we calculate the numerical solution of (|l.ip . (II. 2p 
for $(u) = — (1 — w)^ and the degenerating diffusion coefficient 



Aiu) 







for < u < 10, 



O.l(u-lO) foru>10. 
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Ax 




conv. 
rate 




conv. 
rate 




conv. 
rate 




conv. 

rate 


0.020 


0.239 


- 


0.317 


- 


0.915 


- 


0.695 


- 


0.010 


0.133 


0.845 


0.146 


1.122 


0.513 


0.834 


0.442 


0.655 


0.005 


0.061 


1.135 


0.069 


1.070 


0.246 


1.062 


0.200 


1.144 


0.004 


0.048 


1.018 


0.054 


1.090 


0.181 


1.369 


0.164 


0.891 


0.002 


0.021 


1.168 


0.024 


1.161 


0.082 


1.150 


0.073 


1.163 


0.001 


0.008 


1.360 


0.009 


1.399 


0.036 


1.167 


0.032 


1.200 



Table 1. Example 1: Numerical error at ti = 0.1 and ^2 = 0.25. 

The initial datum is given by 

[5 for 0.1 < x < 0.2, 7 for 0.8 < a; < 0.9, 
uo(x) — < 

[8 for 0.6 < x < 0.7, otherwise. 

Note that Co = 2 in Example 1, where Co is defined in (jl.ll|) . and v* — CqJ2 — 1, 
so that the function $ corresponds to (|1.12|) . where the constant of integration 
is —1, and that uq is chosen such that (|1.8I) is satisfied. Moreover, in our case 
^"{v) = -2 < 0, and $(0) = $(Co) = -1. Nagai and Mimura [13] show that 
under these conditions, and for the integrated diffusion coefficient given by (|1.16p . 
the solution converges in time to a compactly supported, stationary travelling-wave 
solution, which represents the aggregated group of individuals and is defined by the 
time- independent version of (|l.ip . 

In Figure [1] we show the numerical approximations for v and u for < t < 0.5 
and for Ax = 0.001. As predicted, for each fixed time the data {«"} are monoton- 
ically increasing, and the numerical solution for u indeed displays the aggregation 
phenomenon, and terminates in a stationary profile, even though the assumptions 
on ^{•) stated in [53] are not satisfied here. This supports the conjecture that 
a similar travelling wave analysis can also be conducted in the present strongly 
degenerate case, to which we will come back in a separate paper. 

In Tabled] we show the error at t\ =0.1 and t^ — 0.25 in the L^ norm for u 
(denoted as e^', i — 1,2) and in the L°° norm for v (denoted as e^*, i = 1,2), 
where we take as a reference the solution calculated with Aa; — 0.0002. We find an 
experimental rate of convergence in both cases greater than one. For small Ax this 
behavior is possibly related to the proximity of the reference solution. One should 
expect a real order of convergence at most one since the w-scheme is monotone. 

In Figure [5] we compare the numerical approximations for v and a for different 
mesh sizes at the simulated time t — 0.25. 

5.2. Example 2. This example represents a slight modification of Example 1, 
namely we choose j4(-) and uq as in Example 1 but use 

\-(\~vf forO< w< 1, 
^ ' [-(l-v)2 forw> 1. 

The function $ has its maximum in u* = 1 and satisfies $(0) = "I>(Co), as in 
Example 1, so we expect to see an aggregation phenomenon and the formation of 
a stationary travelling wave, even though here $ is not symmetric with respect to 
V* . In Figure [3] we show the numerical approximation of u for < t < 0.5 and 
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Figure 2. Example 1: Numerical approxim.ation of w (top) and u 
(bottom) for several mesh sizes at t = 0.25. 



0.75 



Ax — 0.001. The solution behavior is similar to that of Example 1, but the shapes 
of the "peaks" are slightly different, in particular the final "peak" is asymmetric. 
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Figure 3. Example 2: Numerical approximation of u, obtained 
via (SSI with Ax = 0.001. 



Ax 




conv. 
rate 




conv. 
rate 




conv. 
rate 




conv. 
rate 


0.020 


0.167 


- 


0.330 


- 


0.329 


- 


0.297 


- 


0.010 


0.083 


1.010 


0.166 


0.994 


0.185 


0.834 


0.161 


0.884 


0.005 


0.039 


1.099 


0.079 


1.066 


0.105 


0.812 


0.086 


0.912 


0.004 


0.031 


0.932 


0.062 


1.097 


0.089 


0.733 


0.064 


1.322 


0.002 


0.014 


1.195 


0.028 


1.170 


0.043 


1.060 


0.034 


0.886 


0.001 


0.006 


1.165 


0.010 


1.414 


0.021 


1.056 


0.016 


1.115 



Table 2. Example 3: Numerical error for u and v aX ti =0.1 and 
t2 = 0.25. 



5.3. Example 3. We now choose <I> and uq as in Example 1, but define A{-) by 



A{u) 



0.05u 


for < M < 5, 


0.25 


for 5 < M < 10, 


0.05m- 


~ 0.25 for u > 10. 



Figure [5] shows the results for Ax = 0.001 and t G [0,0.5]. Again, a stationary 
single-peak solution is forming, including a jump between m = 5 and u — 10, in 
agreement with the flatness of A{u) for u S [5, 10]. Table [2] shows the error of the 
approximations of v and u at ti = 0.1 and ^2 = 0.25. The reference solution has 
been calculated with Ax = 0.0002. The observed rate of convergence is again one. 
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Figure 4. Example 3: Numerical approximation of u, obtained 
via ESI) for Ax = 0.001. 



5.4. Example 4. In this example we utilize a flux function with several extrema 
given by $(w) = — 0.5(cos(w7r) + l) combined with the integrated diffusion coefficient 



Aiu) 



f for < M < 10, 

[O.l(u-lO) foru>10 



and the initial datum 



10 for xe [0.05,0.15], 9 for a; e [0.6, 0.7], 
uo{x) == { 14 for X £ [0.3,0.5], 8 for a; e [0.9, 1], 

otherwise. 



The result is shown in Figure [5] for Ax = 0.001. We observe the formation of 
three groups, but the third moves to the right "looking for more" mass since it is not 
a full state, in the sense of the Nagai and Mimura [23 condition for the formation of 
stationary travelling waves. In addition to Figure [S] we show in Figure [5] a contour 
plot of the numerical approximation of v for this example. The contour lines of v 
correspond to trajectories of "individuals" . Table |3] shows the error for v and u 
taking as a reference the solution calculated with Ax = 0.0002. We again find the 
order of convergence predicted for monotone schemes. 
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Figure 5. Example 4: Numerical approximation of u, obtained 
via K^ for Aa; = 0.001. 
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Figure 6. Example 4: Contour lines of the numerical approxima- 
tion of « for Aa; = 0.001. 
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Ax 




conv. 
rate 


e*2 


conv. 
rate 




conv. 
rate 


e*2 


conv. 
rate 


0.020 


0.361 


- 


0.381 


- 


1.420 


- 


1.244 


- 


0.010 


0.189 


0.933 


0.201 


0.923 


0.892 


0.671 


0.709 


0.811 


0.005 


0.095 


0.992 


0.101 


0.994 


0.509 


0.809 


0.356 


0.993 


0.004 


0.080 


0.771 


0.080 


1.048 


0.398 


1.101 


0.262 


1.374 


0.002 


0.042 


0.939 


0.040 


0.981 


0.216 


0.883 


0.145 


0.857 


0.001 


0.019 


1.130 


0.020 


1.000 


0.104 


1.047 


0.072 


1.000 



Table 3. Example 4: Numerical error for u and v at ti =0.1 and 
i2 = 0.25. 
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Figure 7. Example 5: Numerical approximation of u, obtained 
via (HSl) for Aa; = 0.001. 



5.5. Example 5. Here we calculate the numerical approximation of u for A{-) as 
in Example 4, but with $ and uq given by the respective equations 



$(w) 



uo{x) 



I -0.5(cos(i;7r) + 1) for < w < 2, 
[(w- 2)2-1 for w > 2, 

14 for X e [0.15, 0.3], 18 for x e [0.8, 0.95] 
17 for xe [0.6,0.7], otherwise. 



In Figure[7]we show the result for Ax — 0.001. We see that the spare mass (i.e. 
the mass that can not get in the first group) "dilutes" to the right. 

5.6. Example 6. We consider now the same initial data and parabolic term A{-) 
as in Example 1, but employ a function $ with several extrema given by $(w) — 
—0.5 (cos(27ri') + 1). Accordingly with the results of Example 1 we expect a steady 
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Figure 8. Example 6: Numerical approximation of u, obtained 
via ESI) for Ax = 0.001. 



state consisting of two traveling waves since $(0) = <J>(Co). Figure [5] shows the 
numerical result for u for Ax = 0.001, which confirm our claim. 
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